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We review recent developments in the understanding of meson properties as solutions of the Bethe-Salpeter 
equation in rainbow-ladder truncation. Included are recent results for the pseudoscalar and vector meson masses 
and leptonic decay constants, ranging from pions up to cc bound states; extrapolation to bb states is explored. 
We also present a new and improved calculation of F 7r {Q 2 ) and an analysis of the ^77 transition form factor for 
both 7r(140) and 7r(1330). Lattice-QCD results for propagators and the quark-gluon vertex are analyzed, and the 
effects of quark-gluon vertex dressing and the three-gluon coupling upon meson masses are considered. 
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1. DYSON SCHWINGER 
OF QCD 



EQUATIONS 



The Dyson-Schwinger equations [DSEs] are the 
equations of motion of a quantum field theory. 
They form an infinite hierarchy of coupled inte- 
gral equations for the Green's functions (n-point 
functions) of the theory. Bound states (mesons, 
baryons) appear as poles in the Green's functions. 
Thus, a study of the poles in n-point functions us- 
ing the set of DSEs will tell us something about 
hadrons. For recent reviews on the DSEs and 
their use in hadron physics, see Refs. |1I2I3| . 

1.1. Quark propagator 

The exact DSE for the quark propagator is 1 

S{p)~ l —i i>Z 2 + m q {p) Z4 

+ Zi f g 2 D^(q)^S(k)Ti(k,p), (1) 

Jk 

where D^ v (q = k — p) is the renormalized dressed 
gluon propagator, and T\,(k,p) is the renormal- 
ized dressed quark-gluon vertex. The notation J, 

stands for J A d 4 fc/(27r) 4 . For divergent integrals 
a translationally invariant regularization is neces- 
sary; the regularization scale A is to be removed 
at the end of all calculations, after renormaliza- 
tion, and will be suppressed henceforth. 



x We use Euclidean metric {7 M ,7i/} 
a ■ b = aibi - 
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7m = 7m and 



S( P ) 



The solution of Eq. can be written as 
Z{P 2 ) 



i i> + M (p 2 



(2) 



renormalized according to S(p) = ij/> + rn(^) 
at a sufficiently large spacelike /i 2 , with m(p) the 
current quark mass at the scale fi. Both the prop- 
agator, S(p), and the vertex, depend on the 
quark flavor, although we have not indicated this 
explicitly. The renormalization constants Z2 and 
Z4 depend on the renormalization point and on 
the regularization mass-scale, but not on flavor: 
in our analysis we employ a flavor-independent 
renormalization scheme. 

1.2. Mesons 

Bound states correspond to poles in n-point 
functions: for example a meson appears as a 
pole in the 2-quark, 2-antiquark Green's function 
G"( 4 ) — (0|<7i<72<fi<Z2|0). In the vicinity of a meson, 
i.e. in the neig hborhood of P 2 = -M 2 with M 
being the meson mass, such a Green's function 
behaves like 



G (4) 



x(p ,Pt;P)x(h,k ;P) 
P 2 + M 2 



(3) 



where P is the total 4-momentum of the meson, 
p and pi are the momenta of the outgoing quark 
and incoming quark respectively, and similarly for 
ki and k . Momentum conservation relates these 
momenta: p — pi = P = k Q — ki. 
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The function xiPoiPi\ P) describes the projec- 
tion of the bound state onto its qq Fock-space 
component 2 (sometimes referred to as the Bethe- 
Salpeter [BS] wavefunction, which is different 
from the light-cone wave function). It satisfies 
the homogeneous Bethe-Salpeter equation [BSE] 

T(p ,pi;P)= K(p ,pi;k ,ki)x(ko,ki;P) , (4) 

Jk 

with r the Bethe-Salpeter amplitude [BSA] 

T(k , h-P) = Siko)' 1 x(k , h-P) Sihr 1 . (5) 

The kernel K is the qq irreducible quark- 
antiquark scattering kernel. This homogeneous 
integral equation has solutions at discrete values 
P 2 = —M 2 of the total meson 4- momentum P. 

In order to identify the meson mass M from the 
BSE, Eq. (@J, it is convenient to first introduce 
the linear eigenvalue A(P 2 ) of the kernel 



T(p ,pf,P) 
X(P 2 



K(p 0l pi\k ,ki) x(k ,h;P) , (6) 



which has solutions for all values of P 2 . Since we 
use a Euclidean metric, we have to analytically 
continue this eigenvalue equation to the timelike 
region, i.e. negative values of P 2 , and then find 
M such that A(— M 2 ) = 1. The corresponding 
eigenvector T(p ,pi; P) is then a solution of the 
original BSE, Eq. QJ, i.e. it is the BSA of a 
qq bound state. This meson BSA is normalized 
according to 



2P„ 



Nr. 



d 



k , q 



Tr[x(fci, k ) K(k~ , ki;q Q , q t ) x(q 0> <Zi)] 



Tr[f (ki, k Q ) S(k ) T{k , h) S(h)] 



Q=P 

(7) 



at P 2 — —M 2 , with k D — ki = P = q a — qi and 
k Q k Q — Q — q Q q~i. 

The properly normalized BSA T(p ,pi] P) [or 
equivalently, x(Po,Pf, P)} completely describes 

2 Here and subsequently, q refers to dressed quark (quasi- 
particle) states. 



the meson as a qq bound state. Different types of 
mesons, such as pseudoscalar or vector mesons, 
are characterized by different Dirac structures. 
The ground state in any particular spin-flavor 
channel corresponds to the largest eigenvalue A 
in that channel being one. The first excited state 
is determined by the next-largest eigenvalue A be- 
ing one, and so on for higher excited states. 

2. RAINBOW-LADDER TRUNCATION 

A viable truncation of the infinite set of 
DSEs has to respect relevant (global) symme- 
tries of QCD such as chiral symmetry, Lorentz 
invariance, and renormalization group invariance. 
For electromagnetic interactions we also need 
to respect current conservation. The so-called 
rainbow-ladder truncation respects these proper- 
ties. In this truncation, the kernel of the meson 
BSE is replaced by an (effective) one-gluon ex- 
change 



K(Poj Pi] k , ki) 



-Ai: a {q 2 )D^\q)^®^ lv . 



(8) 



where q—p — k = pi — ki, and a(q 2 ) is an effec- 
tive running coupling. The corresponding trun- 
cation of the quark DSE is 



Z ig 2 D^(q)ri(k lP )^ 

A^a{ q 2 )D^{q) lv ^. 



(9) 



This truncation is the first term in a system- 
atic expansion 0j of the quark-antiquark scat- 
tering kernel K; asymptotically, it reduces to 
leading-order perturbation theory. Furthermore, 
these two truncations are mutually consistent in 
the sense that the combination produces vector 
and axial-vector vertices satisfying their respec- 
tive Ward identities. In the axial case, this en- 
sures that in the chiral limit the ground state 
pseudoscalar mesons are the massless Goldstone 
bosons associated with chiral symmetry break- 
ing 5 6. In the vector case, this ensures, in 
combination with impulse approximation, elec- 
tromagnetic current conservation [Jj. We will 
come back to this point later, when we discuss 
electromagnetic processes. 
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Figure 1. The effective ladder-rainbow kernel 
from the Maris-Tandy [MT] model |5] compared 
to LO perturbative QCD and to that deduced 
from quenched lattice QCD data on the gluon 
and quark propagators [Jj]. 



Figure 2. Dynamical quark mass function using 
the rainbow-ladder truncation of Ref . 151 ■ 



the one-loop perturbative running quark mass 



M(p 2 ) 



The ultraviolet behavior of the effective run- 
ning coupling is dictated by the one-loop renor- 
malization group equation; the infrared behav- 
ior of the effective interaction is modeled, and 
constrained by phenomenology. Here, we are us- 
ing a 2-parametcr model for the effective inter- 
action |SJ, fitted to give a value for the chiral 
condensate of (240 MeV) 3 and / w = 131 MeV. 
This effective running coupling, with Aqcd = 
0.234 GeV and Nf = 4, is shown in Fig. □ and 
agrees perfectly with pQCD for q 2 > 25 GeV 2 . 
The model is finite in the infrared region, and 
it is very similar to an effective interaction that 
was recently deduced [5] from quenched lattice 
QCD results for the gluon and quark propaga- 
tors. Indications are that unquenching does not 
change the infrared behavior of the running cou- 
pling much ^U], and thus we believe that this 2- 
parameter model is a realistic parametrization of 
the effective quark-quark interaction in the space- 
like region. 

2.1. Quark propagator 

A momentum-dependent quark mass function 
M(p 2 ) is an essential property of QCD. In the 
perturbative region this mass function reduces to 



In 



p 2 /A 2 



QCD 



(10) 



with 7 ?n = 12/(33 — 2Nf) the anomalous mass 
dimension. Dynamical chiral symmetry breaking 
[D%SB] means that this mass function is nonzero 
even when the current quark mass is zero. In the 
chiral limit the mass function is |llj 



M chilal (p 2 ) 

27T 2 7 m 



W) 



iln 



p 2 /A 2 



QCD 



1-7* 



(11) 



with (qq) the renormalization-point-independent 
vacuum quark condensate [Hi- 
lt is a longstanding prediction of DSE studies 
in QCD that both the mass function M{p 2 ) and 
the wavefunction renormalisation function Z(p 2 ) 
receive strong momentum-dependent corrections 
at infrared momenta, see e.g. Refs. |1I2I3| and 
references therein. Provided that the (effective) 
quark-quark interaction reduces to the perturba- 
tive running coupling in the ultraviolet region, it 
is also straightforward to reproduce the asymp- 
totic behavior of Eqs. and JITJ |12ll3j . Both 
these phenomena are illustrated in Fig. [21 using 
the model for the interaction of Fig. Here we 
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can also see that in the infrared region the dynam- 
ical mass function of the u and d quarks becomes 
very similar to that of chiral quarks. This is a di- 
rect consequence of DxSB, and leads to a dynam- 
ical mass function of several hundred MeV for the 
light quarks in the infrared region, even though 
the corresponding current quark masses are only 
a few MeV. Thus we connect the current quark 
mass of pcrturbative QCD with a "constituent- 
like" quark mass at low energies. 

These predictions for the quark mass func- 
tion have been confirmed in lattice simulations 
of QCD [14115116] . Pointwise agreement for a 
range of quark masses requires this interaction 
to be flavor-dependent 0, suggesting that dress- 
ing the quark-gluon vertex T1(q,p) is important. 
The consequences of a dressed vertex for the me- 
son BSEs are currently being explored |17| . This 
will be discussed in more detail in Sec. El 

2.2. Pseudoscalar and vector mesons 

The most general structure of a pseudoscalar 
meson BSA is 

T PS (k + T 1 P,k-(l-r 1 )P) = 

j 5 [iE+ PF+ fbG + a^kpPuH], (12) 

where the invariant amplitudes E, F, G and H 
are Lorentz scalar functions of k 2 and k ■ P, and 
they also depend on the choice for rj. The natural 
choice for mesons with equal-mass constituents 
(like a pion) is rj — i , though physical observables 
are of course independent of 77. 

In the chiral limit, the axial- vector Ward- 
Takahashi identity [WTI] relates the pseudoscalar 
meson BSA to the (inverse) quark propagator 5^. 
As a consequence, if there is DxSB, there is 
a massless pseudoscalar bound state (Goldstone 
boson) and vice versa: a massless pseudoscalar 
bound state indicates that chiral symmetry is bro- 
ken dynamically. Furthermore, we find 



E{p\p 2 ) 



M(p 2 ) 



(13) 



U z{p 2 ) ■ 

Thus, in the chiral limit we do not have to solve 
any bound state equation in order to obtain the 
coefficient function of the canonical Dirac struc- 
ture of a pseudoscalar meson. This also shows 
why constituent quark models are inadequate in 
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Figure 3. Pion and excited pion mass and decay 
constant, as function of the current quark mass 
(normalized by the realistic u/d quark masses); 
for comparison we have also included the vector 
meson mass. 



describing pions: for a realistic description of pi- 
ons as qq bound states one necessarily has to have 
a (strongly) momentum dependent quark mass 
function. 

Away from the chiral limit, the axial- vector 
WTI relates the pseudoscalar and pseudovector 
spin projections of the BS wavefunction 



/ps ^ps 
where 



fpsPu 



2m,M rps(n) 



(14) 



^ps(m) = -iZ^Nc / Tr[x(fc 0) fc i )75 



Z 2 N C Tr[ X (k 0l h) l5lfl ], (15) 



(16) 



This is an exact relation in QCD, and holds in- 
dependent of the current quark mass. The Gell- 
Mann-Oakes-Renner relation is a direct corollary 
of this relation, but that is not the only observable 
consequence. It also means that the decay con- 
stant /ps of excited pions vanishes in the chiral 
limit 18 . This relation is indeed satisfied in the 
rainbow-ladder truncation, as can be seen from 
Fig. |3J both for the ground state pion and for 
its first radial excitation Furthermore, for 



5 



Table 1 

DSE results 8 for the meson masses and decay 
constant, together with experimental data |2U] . 





experiment 


calculated 




(estimates) 


(t fitted) 


rn u=d 
''V=lGeV 


4 - 11 MeV 


5.5 MeV 


TO ^=lGcV 


100 - 300 MeV 


125 MeV 


- mi 


(0.24 GeV) 3 


(0.241+) 3 




0.138 GeV 


0.138+ 


u 


0.131 GeV 


0.131+ 


m K 


0.496 GeV 


0.497+ 


Ik 


0.160 GeV 


0.155 


nip 


0.770 GeV 


0.742 


f P 


0.216 GeV 


0.207 


m K * 


0.892 GeV 


0.936 


fx* 


0.225 GeV 


0.241 


m<f, 


1.020 GeV 


1.072 


f& 


0.236 GeV 


0.259 



heavy-light pseudoscalar mesons it predicts ^Hj 
that the decay constant decreases like 1 / \J~M. 

The next-lightest mesons are the vector 
mesons. The BSA of a massive vector meson can 
be decomposed into eight Dirac structures 

r^(fc + v p,k - (l - v)P) = 

8 

^r(fc 2 ,fc-F;r,)T;(fc,P), (17) 

where T^(k,P) are eight independent transverse 
Dirac tensors. Again, the invariant amplitudes /* 
are Lorentz scalar functions of k 2 and k ■ P, and 
depend on r\. 

With the model of Ref. [B], we obtain good 
agreement with the experimental values for the 
light pseudoscalar and vector meson masses and 
leptonic decay constants, see Table ^ These 
results show little sensitivity to variations in 
the model parameters, as long as the integrated 
strength of the effective interaction is strong 
enough to generate an acceptable amount of chi- 
ral symmetry breaking, as indicated by the chiral 
condensate. This is not true for heavier states: 



e.g. the radially excited pion is quite sensitive to 
details of the interaction [2] • 

2.3. Frame independence 

The BSE is usually solved in the rest-frame of 
the meson. To be explicit, the most convenient 
frame is that characterized by P M = (i M, 0, 0, 0). 
However, this is not required: One of the advan- 
tages of the DSE approach to hadron physics is 
its explicit Poincare covariance. This makes it a 
particularly useful tool for studying electromag- 
netic form factors and other processes. No mat- 
ter what frame one chooses for calculating say the 
pion electromagnetic form factor, at least one of 
the pions is moving. 

As a test, we have calculated the static w and 
p properties in a moving frame P M — (i E, q, 0, 0) 
where q is the 3-momentum of the moving meson. 
Within this frame we solve again the eigenvalue 
equation, Eq. |J?J|, search for P 2 = — M 2 such 
that A(P 2 ) = 1, and calculate the corresponding 
electroweak decay constant. Numerically, this is 
a bit of a tour-de-force, since the Lorentz scalar 
functions of Eqs. I|12(l and 1)17(1 are now functions 
of a radial variable k 2 and two angles 

k-P — i k E cos a + k q sin a cos/3 , (18) 

and the integral equation has to be solved in 
the three independent variables k 2 , a, and (3. 
With current computer resources, this can be 
done without further approximations, and the re- 
sults, shown in Fig. 0] arc indeed independent of 
the meson 3-momentum, illustrating that this ap- 
proach is indeed Poincare covariant. We can now 
use this same approach to calculate meson form 
factors in an explicitly covariant manner. 

3. ELECTROMAGNETIC PROCESSES 

In order to study electromagnetic processes, we 
need to know the dressed quark-photon coupling, 
in addition to the meson BSAs and the dressed 
quark propagators. 

3.1. Quark-photon vertex 

The g<77 vertex is the solution of the renormal- 
ized inhomogeneous BSE with the same kernel K 
as the homogeneous BSE for meson bound states. 
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Figure 4. Pion and p mass and decay constant 
calculated in a moving frame, as function of the 
meson 3-momentum. 



That is, for photon momentum Q, we have 



r M (p ,Pi) 



^2 7m+ / K(p ,pi;k 0l ki) 

Jk 

xSi^r^kJSiki), (19) 



with p and pi the outgoing and incoming quark 
momenta, respectively, and similarly for k and 
ki, with p — pi = k a — fcj = Q. Because of gauge 
invariance, the vertex satisfies the vector WTI 



) = S- 1 (p + )-S-\p.). 



(20) 



The most general form of the quark-photon 
vertex T^ipcPf, Q) requires a decomposition into 
twelve Dirac structures. Four of these covariants 
represent the longitudinal components which are 
completely specified by the WTI in terms of the 
(inverse) quark propagator and they do not con- 
tribute to elastic form factors. The transverse 
vertex can be decomposed into eight Dirac ten- 
sors T^(p; Q) with the corresponding amplitudes 
being Lorentz scalar functions. 

Note that solutions of the homogeneous version 
of Eq. i|T§|) define vector meson bound states with 
masses M v = — Q 2 at discrete timelike momenta 
Q 2 . It follows that T M has poles at those locations 
and, in the neighborhood of Q 2 = — My, behaves 



like 



r + M 2 



(21) 



where is the vector meson BSA, and /y the 
electroweak decay constant. The fact that the 
dressed qqj vertex exhibits these vector meson 
poles lies behind the success of naive vector- 
meson-dominance [VMD] models; the effects of 
intermediate vector meson states on electromag- 
netic processes can be unambiguously incorpo- 
rated by using the properly dressed qqj vertex 
rather than the bare vertex 7 M [21] . 

3.2. Impulse approximation 

The generalized impulse approximation allows 
electromagnetic processes to be described in 
terms of dressed quark propagators, bound state 
BSAs, and the dressed grjy-vertex. Consider for 
example the 3-point function describing the cou- 
pling of a photon with momentum Q to a pseu- 
doscalar meson ab, with initial and final momenta 
P ± Q/2, which can be written as the sum of two 
terms 



A^(P,Q) = Q a A° ab + Q h A™ 



(22) 



where Q is the quark or antiquark electric charge, 
and where A ala {P, Q) and A ahh {P, Q) describe the 
coupling of a photon to the quark (a) and anti- 
quark (6) respectively. In impulse approximation, 
these couplings are given by 

Kf{P,Q) = iN c jTr[Tl(q_, q+ )x al (q + ,q) 

x^Cgr 1 * 5 "^-)] , (23) 

with q=k- P/2 and q± = k + P/2 ± Q/2, and 
similarly for A" 66 . The corresponding meson elas- 
tic form factor is defined by 



Af(P,Q) = 2P^F(Q 2 ) 



(24) 



3.3. Current conservation 

Electromagnetic current conservation dictates 
F(0) = Q a + Q b and Q^AfiPQ) = 0. It is 
straightforward to show that both these condi- 
tions are satisfied in impulse approximation, in 
combination with rainbow-ladder truncation. At 
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Q = the q<77 vertex is completely specified by 
the differential Ward identity 



d 

*r M (p,p) = - — s~ (p) 



(25) 



If this is inserted in Eq. Ij23(l . one finds after a 
change of integration variables 



2P fl F 



aab 



(0) 



N r /Tr 



T ab (q',q)S b (q)T ba (q, 



dS a (q') 

8Pa 



(26) 

with q' = q + P. Comparing this expression 
to Eq. l(7|l. we recognize that the requirement 
F aafc (0) = 1 follows directly fr om the canonical 
normalization condition for r ab provided that (1) 
the BSE kernel K is independent of the total me- 
son momentum P and (2) the aaj vertex satisfies 
the differential Ward identity [7]. 

The second constraint, Q^A^ b (P,Q) = 0, is 
satisfied if the quark-photon vertex satisfies the 
vector WTI. Inserting Eq. JSHJl into Eq. J23|), we 
find 



Q 



aab ( 



(P,Q)=N C Tr[ X ah (q + , q)f ba (q,q^] 



N, 



[TT[T ah ( q+ ,q)x ha (q,q-)} , 
Jk 



(27) 



which clearly vanishes. 

The ladder BSE kernel is indeed independent of 
the meson momentum P, and with quark propa- 
gators dressed in rainbow approximation and the 
dressed quark-photon vertex calculated in ladder 
approximation, the both the vector WTI and the 
differential Ward identity are satisfied [221; thus 
impulse approximation for the quark-photon ver- 
tex in combination with rainbow-ladder trunca- 
tion is consistent in the sense that the correct 
electric charge of the meson is produced indepen- 
dent of the model parameters and the resulting 
meson electromagnetic current is conserved. 

3.4. Pion form factor 

For the pion electromagnetic form factor, work- 
ing in the isospin limit, we have 

2P M iv(Q 2 ) - iN c [Tr[ X n (q,q+)T u (q + ,q_) 

Jk 
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Figure 5. Our results for Q 2 Ft T (Q 2 ), with es- 
timates of our numerical error bars, compared 
to a VMD model and to experimental data from 
CERN ID and JLab |25j. 



with the propagators solutions of the quark DSE 
in rainbow truncation. The pion BSAs and the 
qepi vertex are solutions of their respective BSEs 
in ladder truncation, using exactly the same mo- 
mentum frame as used in the form factor calcu- 
lation. That does mean re-calculating the pion 
BSAs for every value of the photon Q 2 , but the 
advantage is that one does not have to do any 
interpolation or extrapolation on the numerical 
solutions of the BSE. Our results for spacelike Q 2 
are shown in Fig. |SJ They are in excellent agree- 
ment with the available experimental data ,24 25 
and with our previous results 3 . 

Technical difficulties due to the appearance 
of a pair of complex conjugate singularities in 
the quark propagator limit our calculation to 
Q 2 < 5 GeV 2 . Up to this value of Q 2 we see 
no significant deviation from a naive monopole 
behavior in our calculation of the form factor. 
It would be very interesting to see whether or 
not TJNAF can establish a clear deviation from 
a monopole behavior, because pQCD predicts a 



Xx«(q_,q)S(q)- 



(28) 



3 The current calculations improve on Refs. |22l2.'il by elim- 
inating: the Chebyshev expansion for the angular depen- 
dence of the restframe BSE solutions, and the interpola- 
tion and extrapolation techniques used to connect frames. 



8 



• — • DSE calculation 
x * lattice Bonnet 
+■■+ lattice v.d. Heide 
* expt. 

-- VMD: 6/M 2 



^ 0.3 - 



0.2 0.4 0.6 

m„ [GeV] 



Figure 6. Pion charge radius as function of the 
quark mass within our DSE model, and as ex- 
tracted from lattice simulation |2fil27| . 



significantly lower value for Q 2 Ft T (Q 2 ) at (asymp- 
totically) large Q 2 . 

3.5. Quark mass dependence 

We have calculated Fk(Q 2 ) for a range of cur- 
rent quark masses. For all masses considered, 
the form factor behaves like a monopole, with a 
monopole mass that is slightly smaller than the 
p-meson mass. A similar behavior has been found 
in quenched lattice simulations 26 27 . 

The obtained charge radius, defined by 
r 2 = —6F'(0), is also in excellent agreement with 
the experimental value. Our calculated charge ra- 
dius as a function of the quark mass is shown in 
Fig. El together with the charge radius obtained 
from lattice simulations. For comparison, we also 
show the charge radius one would obtain from a 
naive VMD model, using the corresponding (cal- 
culated) vector meson mass. Again, both our nu- 
merical result for r 2 and that of lattice simula- 
tions are represented reasonably well by a VMD 
model for the charge radius. 

An effect that is missing both from our calcu- 
lations and from quenched lattice simulations are 
pion loop corrections. At present it is not clear 
how to incorporate meson loops self-consistently 
in both the DSE/BSEs and in the approxima- 
tion for the photon-hadron coupling. However, 



in Ref. it was demonstrated that the dressed 
quark core generates most of the pion charge ra- 
dius, and that pion loops contribute less than 15% 
to r 2 at the physical value of the pion mass. For 
larger values of the current quark mass (and thus 
of the pion mass), we expect these corrections to 
become negligible. Also for larger spacelike val- 
ues of Q 2 the effects from meson loops decrease, 
and for Q 2 > 1 GeV 2 we expect the contribution 
of such loops to be negligible. 

3.6. 7T7 7 transition form factor 

Using impulse approximation, the coupling of 
a pion to two photons is, in the isospin limit 



A^(Qi,Q 2 ) 



Ti[x*(q,p)rl(p,k)S(k)TZ(k,q)] , (29) 



k 

2 i a e 

7T 



/3 QlaQ2pGit~ 11 {Q\, Q\) , (30) 



where the internal momenta p, q are determined 
in terms of the integration momentum k and the 
external momenta Qi,Q2 by momentum conser- 
vation; furthermore, the pion is on-shell, but the 
photons not necessarily. 

For on-shell photons, Q 2 = = Q 2 ,, the partial 
decay width ir° — > 77 is given by 



r 

* 7T 



77 



2 3 

16tt 3 



Gwy(0,0) 



(31) 



This decay width is governed by the axial 
anomaly as prescribed by electromagnetic current 
conservations and chiral symmetry. In the chiral 
limit, the axial anomaly gives 

G 7r77 (0, 0) = — « 5.4 GeV" 1 , (32) 

2 fir 

and the resulting decay width 



7T77 



2 3 

a cm m n 



64tt3/ 2 



(33) 



is within 2% of the experimental width of 7.8 eV. 
The corrections due to finite pion mass are small. 

4 This expression is conventionally given in terms of 
fn = fir/v^t the pion decay constant in the convention 
where its value is 92 MeV, rather than the convention 
/tt = 131 MeV used throughout this work. 
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Figure 7. Symmetric ^77 transition form factor 
G(Q 2 , Q 2 ) for both the ground state pion and its 
first radial excitation. Adapted from \21\ . 



Figure 8. Asymptotic behavior of the symmetric 
7T77 transition form factor. Adapted from 21 . 



The axial anomaly is preserved by the ladder- 
rainbow truncation of the DSEs combined with 
an impulse approximation for the 7r°77 vertex 
because the relevant manifestations of electro- 
magnetic gauge invariance and chiral symmetry 
are present |7I29| . We do indeed reproduce the 
above value for this form factor, as can be seen 
from Fig. [7| Our results also agree with the 
experimental data for the transition form factor 

G ff7 * 7 (<9 2 ,o) H3E2- 

On the other hand, the radially excited pion 
effectively decouples from the axial anomaly |2T] . 
and its coupling to two on-shell photons is not 
constrained by the anomaly. A naive application 
of Eq. (pf2"|l to the first radially excited pion would 
suggest an enormously large value of 625 GeV -1 
for its coupling G(0, 0) to two photons, because its 
decay constant is very small, w —0.0016 GeV, 
and in the chiral limit this decay constant van- 
ishes, as is evident from Eq. l(T3|l and Fig.|31 How- 
ever, numerically we find that for the excited pion 
this coupling is about 0.7 GeV -1 . Furthermore, 
the chiral limit 71-77 form factor is at low-Q 2 al- 
most identical to the curves shown in Fig. [7| both 
for the ground state and for the excited pion. 

The asymptotic behavior of this form factor is 
given by the lightcone operator product expan- 



sion [22] 

G{Q\,Ql) -> 27T 2 / 7rn I ^ + q2 

+ °(t'w^w)}-< m » 

where cj = {Q\ — <9l)/( < 5i + Q2) ^ s the photon 
asymmetry. In lightfront QCD J(oj) is related 
to the leading twist pion distribution amplitude 

4> n (x) via 

J(W) = H \-J(2 X -1) ^ X) - (35) 

The normalization of cf) w (x) immediately gives 
J(0) = I for the case of equal photon virtuality. 

Within the present DSE-based approach, it is 
straightforward to reproduce this leading-order 
asymptotic behavior analytically, not only for the 
ground state pion but also for its excited 
states pi] . Also numerically we reproduce 

G(Q 2 ,Q 2 ) - (36) 

as can be seen from Fig. [S] 

The interesting question is what happens with 
the form factor of the (first) radially excited 
pion in the chiral limit. With a nonzero current 



10 



quark mass, this form factor indeed behaves like 
Eq. (|36[1 . and it approaches zero from below as 
Qi — Q\ 00 ' since its decay constant is nega- 
tive. In the chiral limit however, the decay con- 
stant f 7Tl is zero, as dictated by Eq. I|14[) . and 
evident in Fig. |3| It turns out that in that case 
the form factor vanishes like l/Q 4 , with a positive 
coefficient: the form factor has a zero-crossing in 
the spacelike region. The coefficient for this l/Q 4 
behavior depends on the details of the model and 
its parameters, but it is not simply related to f ni . 
The leading order behavior of Eq. however 
is model-independent. 



4. MORE MASSIVE MESONS 

The rainbow-ladder truncation of the set of 
DSEs, using the model interaction of Rcf. 8 , has 
been very successful in describing a wide range of 
light meson properties J2j ■ One question that nat- 
urally arises: How does it do for heavier mesons? 

4.1. Heavy quark mesons 

Let us first address the question of the quark 
mass dependence of the pseudoscalar and vector 
meson properties. In Fig. [5] we show the meson 
mass and decay constants as function of m 9 (19), 
the current quark mass at the renormalization 
point fi = 19 GeV, both for heavy-light mesons 
(light meaning u/d quarks and s quarks) and for 
heavy-heavy mesons. 

From this figure it is apparent that all of the 
meson masses scale approximately linear with the 
current quark mass above about m q ~ 0.5 GeV rs 
6to s , with m s being the strange quark mass. The 
pseudoscalar meson masses show a clear curva- 
ture below this mass scale, as one would expect 
based on the Gell-Mann-Oakes-Renner relation. 
Also the vector meson masses show a slight cur- 
vature at small current quark masses. The re- 
sults for the D, D s , rj c , and J/^> mesons are 
in remarkably good agreement with experiments: 
the calculated masses with m c (19) = 0.88 GeV 
are within 2% of the experimental values. Com- 
pared to the light current quark masses, we find 
Tn c /m u = 238 and m c /m s = 10.5. 

Both the vector and the pseudoscalar masses 
can be fitted reasonably well over the entire range 
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Figure 9. Meson mass as function of m q . Top: 
our calculation (with estimates of numerical er- 
rors) up to m c plus the fit Eq. l|3*7|) : bottom: ex- 
trapolated to mj, using the same fit; the vertical 
dashed lines indicate m u f d , m s , m c , and m^. 



< m q < m c by the same formula 



M l 



Mq + ai(mi + m 2 ) + a 2 (mi + TO2) 2 , (37) 



as indicated by the lines in Fig. Here m, 
are the current quark masses of the constituents 
(at n = 19 GeV); the chiral limit masses are 
Mq S = 0, M y = 0.75 GeV; and the remaining fit 
parameters are af s = 2.96 GeV, = 3.24 GeV, 
and a 2 = 1.12 both for pseudoscalar and vector 
mesons. 

If we use this fit to extrapolate our results to 
the bottom quark mass, we get the bottom panel 
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of Fig. |31 The results for the B, B s , and T mesons 
are in surprisingly agreement with experiments: 
the calculated masses with m/,(19) = 3.8 GeV 
are within 1% of the experimental values. Com- 
pared to the other current quark masses, we find 
mb/m c = 4.3, mb/m s = 45 and mb/m u is about 
1,000. Our prediction for r) b is M Vb = 9.34 GeV 

In Fig. 1 101 shows the corresponding leptonic de- 
cay constants. At small current quark masses, all 
of the decay constants increase with the quark 
mass, but they tend to level off between the s and 
the c quark mass, both for the heavy- light and for 
the heavy-heavy mesons. This is consistent with 
the expected behavior in the heavy-quark limit: 
namely a decrease of the decay constant with in- 
creasing meson mass (or equivalently, increasing 
heavy-quark mass) like / ~ 1/ \/M, at least for 
the heavy-light mesons. Our current calculations 
suggest that the onset for this asymptotic behav- 
ior could be as low as the charm quark mass for 
the qu mesons. However, our results for both 
Jd and fo s are about 20% below their experi- 
mental values, f D = 0.222 ± 0.020 GeV gS| and 
f Dg = 0.266 ± 0.032 GeV ; 2H respectively and 
our result [33] for the decay constant of the J/ip 
is about 25% below its experimental value, ex- 
tracted from the e + e~ decay [201 ■ This could in- 
dicate that the rainbow-ladder truncation is not 
reliable in the charm quark region, and/or our 
model for the effective interaction is not applica- 
ble to charm quarks, despite the fact that the 
meson masses are in good agreement with ex- 
periments. Since the decay constants depend on 
the norm of the BSAs, which in turn depends 
on the derivative of the quark propagators, it is 
not a surprise that the decay constants are more 
sensitive to details of the model than the meson 
masses, and are thus better indicators for defi- 
ciencies in the modeling. 
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Figure 10. Meson decay constant (with estimates 
of our numerical errors) as function of m q , up to 
the charm quark mass region; the vertical dashed 
lines indicate m u /d, rn s , and m c . 



the experimental value |18| . Furthermore, these 
heavier mesons, as well as the scalar mesons |37) . 
tend to be more sensitive to details of the inter- 
action than the p and it. 

This suggests a deficiency of the model and/or 
of the rainbow-ladder truncation for these higher- 
mass states. A possible explanation for this de- 
ficiency lies in the fact that these higher-mass 
states have a different spin structure and (espe- 
cially in case of the radially excited pion) very 
different BS wave functions, and therefore they 
are sensitive to different aspects of the BS kernel 
K. Again, we should study the effects from con- 
tributions beyond the rainbow-ladder truncation 
in more detail in order to resolve these issues. 

5. BEYOND THE RAINBOW 



4.2. Heavier u and d quark mesons 

Another question is that of heavier states con- 
sisting of u and d quarks. Calculations of the a\ 
and b\ mesons using a similar model give results 
for these masses that are several hundred MeV 
too low 



when compared to experiments. 
Also the mass of the first radially excited pion 
is, within this model, significantly smaller than 



Both lattice simulations |38I39| and DSE stud- 
ies [1714014 1] of the quark-gluon vertex indi- 
cate that Tl(q,p) deviates significantly from a 
bare vertex in the nonperturbative region. A 
(flavor-dependent) nonperturbative vertex dress- 
ing could make a significant difference for the so- 
lution of the quark DSE |10ll7j . Simple dressed 
vertex models have indicated that material con- 
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tributions to a number of observables are possi- 
ble with a better understanding of the infrared 
structure of the vertex. These diverse model in- 
dications include: an enhancement in the quark 
condensate |10I17) : an increase of about 300 MeV 
in the b\jh\ axial vector meson mass |36| : and 
about 200 MeV of attraction in the p/u> vector 
meson mass |17) . 

5.1. Lattice-inspired rainbow model 

Over the last five years, the gluon 2-point 
function has become fairly well known through 
explicit lattice-QCD calculations jl^l an d also 
through approximate solutions of the gauge sec- 
tor DSEs |2I1UI43| . A strictly defined rainbow 
truncation of the quark DSE for calculation of 
the dressed quark propagator would proceed from 
this input. However, the modern information on 
the dressed gluon propagator shows a strong in- 
frared suppression and it is not possible to ob- 
tain a realistic value for the condensate (qq) in 
rainbow truncation 03]. In this sense it is clear 
that the empirically successful rainbow-ladder 
kernels developed earlier implicitly include effects 
of quark-gluon vertex dressing for the quark DSE. 
Although such an effective kernel shows a bare 
vertex Dirac matrix structure, the infrared pa- 
rameterization produces a strength and depen- 
dence upon gluon momentum that is over and 
above that of the gluon propagator. 

This point is illustrated in Fig. El 9 , where 
quenched lattice data on the gluon and quark 
propagators are used to deduce a phenomenolog- 
ical quark-gluon vertex by means of the DSE. 
The DSE kernel of Ref. [§| is modeled by 
^a(q 2 )D^(q) -> V(q 2 , m q ) D$(q) in Eq. ©, 
where D l ^(q) represents the quenched lattice re- 
sult for the Landau gauge gluon propagator, and 
V(q 2 ,m q ) is a phenomenological representation 
of vertex dressing whose ultraviolet behavior is 
determined by the requirement that the resulting 
running coupling reproduces the one-loop renor- 
malization group behavior of QCD. Thus, the ver- 
tex dressing has been mapped onto a single ampli- 
tude corresponding to the canonical Dirac matrix 
7^, and depends on the gluon momentum and the 
current quark mass. 

The resulting effective DSE kernel [DSE-Lat] is 



• lattice data 




p [GeV] 



Figure 11. Dynamical quark mass functions ob- 
tained with the DSE-Lat model together with the 
quenched lattice data Adapted from [Sj. 



compared to the MT model in Fig. 2] The param- 
eters are determined by requiring that the DSE 
solutions reproduce the quenched lattice data 
for S(p) in the available domain p 2 < 10 GeV 2 
and m(fi = 2 GeV) < 200 MeV. In this sense, 
the DSE-Lat model represents quenched dynam- 
ics. It is found that the necessary vertex dress- 
ing is a strong but finite enhancement. As one 
would expect, the vertex dressing V(q ,m q ) de- 
creases with increasing m q to represent the effect 
of quark propagators internal to the vertex. 

The model easily reproduces m v with a current 
mass that is within acceptible limits. However the 
resulting chiral condensate (qq) — (0.19 GeV) 3 is 
a factor of 2 smaller than the empirical value 
(0.24 ± 0.01 GeV) 3 . This is attributed to the 
quenched approximation in the lattice data. 

5.2. Modeling the quark-gluon vertex 

Let us denote the dressed quark-gluon vertex 
for gluon momentum k and quark momenta p and 
p+khy ig^-T a (p + k,p). Through 0(g 2 ), i.e., to 
one loop, we have r CT = Zip "f a + r£ + T^ A , see 
Fig. El Here Z\p is the vertex renormalization 
constant to ensure r CT = "f a at the renormaliza- 
tion scale fi. 

The color factor for is (Cp - §Ca) = -JW* 
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Figure 12. The one-loop corrections to the quark- 
gluon vertex dressing. Left: the Abelian-like term 
r A ; Right: the non-Abelian term T^ A . 



whereas the color factor for is \Ck = j • 
This reveals two important issues: The T A con- 
tribution to the quark-gluon vertex is repulsive, 
whereas the non-Abelian term, which in- 

volves the three-gluon vertex, is attractive. Fur- 
thermore, the magnitude of the non-Abelian term 
r^ A is enhanced by a factor N 2 over the T A 
term 5 . 

The term T^ A involves the triplc-gluon vertex 
r^ CT , which is a "new" interaction. However, in 
Landau gauge the corresponding (bare) Slavnov- 
Taylor identity reduces F^ A at zero gluon mo- 
mentum, k = 0, to 

r^ A (p,p) = JlnS (p-q)~i v 

X {^9 2 D (q 2 )^T tw (q)- (38) 

This expression has been used in a recent 
model 0U] f° r the dressed quark-gluon vertex at 
k = 0. It is a non-perturbative extension of the 
two diagrams in Fig. 1121 In the ultraviolet, the qq 
scattering kernel appearing in T A coincides with 

5 Note that the color factor of the color-singlet contri- 
bution to the strong dressing of the quark-photon vertex 

is Cf = N ^ N 1 , i.e. attractive and enhanced by N% over 
the color-octet term: single gluon exchange between a 
quark and antiquark has relatively weak repulsion in the 
color-octet channel, compared to the strong attraction in 
the color-singlet channel. 



the ladder-rainbow kernel; hence the natural ex- 
tension is g 2 D (q 2 ) — > 4ira c ff(q 2 )/q 2 and the ex- 
ternal quark-gluon vertex is taken to be bare. 

For r^ A it takes full advantage of the simpli- 
fication, Eq. (|38|l . that happens for zero momen- 
tum gluons: replace g 2 D (q 2 ) by 4ira e ff(q 2 )/q 2 . 
This term is similar to a contribution to the 
derivative of the quark self-energy but the differ- 
ences are important; at one-loop this term pro- 
vides the explicitly non-Abelian contributions to 
the Slavnov- Taylor identity 05]- The justifica- 
tions for this nonperturbative vertex model are 
consistency and simplicity; no new parameters 
are introduced. 

In Fig. ED we display the results of this model 
in a dimensionless form for comparison with the 
(quenched) lattice data for m(/i) = 60 MeV. The 
general vertex at k = has a representation in 
terms of four invariant amplitudes, but only 3 are 
non-zero here. We write 

r<r(p,p) = i*M(p 2 ) - 4p ff 7 -pX 2 (p 2 ) 

-i2 Pa \ 3 ( P 2 ), (39) 

since the lattice data |38) are provided in terms 
of these Xi(p 2 ) amplitudes. The renormalization 
scale for both the lattice data and the DSE-Lat 
model is fj, = 2 GeV where Ai(/x) = 1. Without 
parameter adjustment, the model reproduces the 
lattice data for Ai and A3 quite well. The lattice 
result for A2, despite the large errors, suggests 
infrared strength that is seriously underestimated 
by the DSE-Lat model. 

The relative contributions to the vertex 
dressing made by T^ A and T A are indi- 
cated by the following amplitude ratios at 
p = 0: Af A /A A = -60, A^ A /A A = -14, and 
A^ A /A A = -12. Thus the non-Abelian term T^ A 
dominates to a greater extent than what the ratio 
of color factors (—9) would suggest. 

Another useful comparison is the correspond- 
ing vertex in an Abelian theory like QED; it is 
given by the differential Ward identity, Eq. (|25|l . 
With S^ 1 (p) = i-y -pA(p 2 ) + B(p 2 ), this leads to 
the correspondance A^ 1 = A, Xf 1 = -A'/2, and 
A^i = B', where /' = df(p 2 )/dp 2 . The Abelian 
Ansatz , while clearly inadequate for the dom- 
inant amplitude Ai below 1.5 GeV, does repro- 
duce the DSE-Lat results for both A2 and A3. 
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Figure 13. The amplitudes of the dressed quark- 
gluon vertex at k = for m q {2 GeV) = 60 MeV. 
Quenched lattice data [HHj are compared to the re- 
sults of the DSE-Lat model jHj, and for \\{p) the 
Abelian Ansatz is also shown. Adapted from 



Figure 14. The relative error of the ladder- 
rainbow truncation relative to the completely re- 
summed model of vertex dressing and correspond- 
ing BSE kernel. Adapted from Ref. [T7|. 



5.3. Effect on meson observables 

There is a known constructive scheme |J that 
defines a diagrammatic expansion of the BSE ker- 
nel corresponding to any diagrammatic expansion 
of the quark self-energy such that the vector and 
axial-vector WTIs are preserved. Among other 
things, this guarantees the Goldstone boson na- 
ture of the flavor non-singlet pseudoscalars inde- 
pendently of model details [SJ . 

This has recently been exploited to produce 
the first study |17| of how an infinite sub-class 
of quark-gluon vertex diagrams, including at- 
traction from three-gluon coupling, contributes 
to vector and pseudoscalar meson masses. The 
symmetry-preserving BSE kernel obtained from 
the self-energy has an infinite sub-class of dia- 
grams beyond ladder; it is manageable here be- 
cause of the algebraic structure afforded from use 
of a modification of the Munczek-Nemirovsky 8- 
function model It is verified that the pseu- 
doscalar qq meson remains a Goldstone boson at 
any order of truncation. As shown in Fig. 1141 the 
ladder-rainbow truncation has at most a 8% re- 
pulsive error as m q is varied up to the b quark 
region. 



Fig. ^] also shows that attraction from the 
dressed vertex is a significant correction to ladder- 
rainbow truncation for the qq light quark vector 
mesons, and reduces steadily in relative impor- 
tance as Tfiq increases. This is consistent with the 
results [§] shown in Fig. ^] where the effective 
vertex dressing V(q 2 ,m q ) required by quenched 
lattice-QCD propagators also decreases with (rel- 
atively light) m q . It is also consistent with our 
results, shown in Fig. [5] for the m q dependence 
of the meson masses from the more realistic MT 
ladder-rainbow model 8 . The reason is that 
the latter model has a phenomcnological infrared 
strength fitted to empirical light quark informa- 
tion, (qq); it therefore implicitly includes much 
of the attraction from vertex dressing. The re- 
sulting m p is only 5% off. The algebraic model 
analysis in Fig. ^] suggests that the relative er- 
ror of a ladder-rainbow vector mass can only get 
smaller with increasing m q ; it should drop by a 
factor of about 5 in going from the u/d quark re- 
gion to the c quark region, and be about 1% for b 
quarks. That is what we find for the MT model 
in Fig. [5| the error is only 2% for c quarks and 
1% for b quarks. Although the algebraic model is 
simple and not all details will be realistic, we ex- 




Figure 15. The two leading-order (1-loop) vertex 
corrections to the rainbow DSE (top) and the cor- 
responding 5 additions to the ladder BSE kernel 
(bottom) to preserve the relevant Ward identities. 
Quark and gluon lines are dressed. 



pect the qualitative aspects of these observations 
to be confirmed whenever a more realistic study 
becomes feasible. 

5.4. Current conservation 

If one goes beyond the rainbow-ladder trunca- 
tion for the propagators, the meson BSAs and 
the quark-photon vertex, one has to go beyond 
impulse approximation for the meson form factor 
to ensure current conservation 23 . For exam- 
ple, following the general procedure of Ref. 0], 
one-loop dressing of the quark-gluon vertex in the 
DSE requires 5 additions to the ladder BSE ker- 
nel to preserve the relevant WTIs, see Fig. IT31 

The resulting BSE kernel K(q 0l g,-; k 0l ki) now 
becomes dependent on the total meson momen- 
tum P = q Q — qi = k a — hi, which means that the 
first term in the normalization condition, Eq.Q, 
is nonzero. With the choice k a — k + P, ki = k, 
this introduces the 4 extra terms in the normaliza- 
tion condition depicted in the top row of Fie . 1161 
These 4 additional diagrams are generated from 
the BSE kernel in the bottom part of Fig. El 
by taking the derivative with respect to the me- 
son momentum P, where P flows through one 
quark propagator only. Since a derivative with re- 
spect to P is equivalent to the insertion of a zero- 
momentum photon according to the differential 
Ward identity, Eq. i|25[l . it is obvious which dia- 
grams have to be added to the impulse approxi- 
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Figure 16. Top: The 4 contributons to the 
BSA normalization condition from the additions 
of Fig. EH to the DSE dynamics. Bottom: The 
corresponding 4 corrections to the impulse ap- 
proximation to the form factor. Slashes denote 
derivatives as described in the text. 



mation to maintain current conservation [23] , and 
they are displayed in the bottom row of Fig. 1161 
The transversality of the photon-meson coupling 
is also maintained by this truncation. 

We note that this procedure systematically ex- 
poses Fock-space components of the meson that 
are included in the exact qq scattering kernel of 
the BSE but which necessarily show up explicitly 
as corrections to the impulse approximation. 

6. CONCLUDING REMARKS 

We have shown how the rainbow-ladder trun- 
cation of the set of DSEs can be applied to study 
pseudoscalar and vector meson properties, from 
pions up to bb bound states. This truncation re- 
spects the relevant symmetries and related Ward 
identities: e.g. the axial-vector WTI, which en- 
sures the Goldstonc nature of the pions in the 
chiral limit; we also preserve the correct QCD 
renormalization properties. Another advantages 
of this approach is its manifest Poincare invari- 
ance, which makes the study of electromagnetic 
form factors relatively simple. In addition to 
the dressed quark propagators and bound state 
amplitudes, one needs the quark-photon vertex, 
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which can be obtained by solving an inhomoge- 
neous BSE [22] • The rainbow-ladder truncation, 
in combination with impulse approximation for 
the form factor, guarantees current conservation. 
Our result for F 7T {Q 2 ) is in remarkable agreement 
with the data from TJNAF - 

The results for the kaon form factors are also 
in good agreement with the data |23| . as are the 
results for the radiative (30831| and strong de- 
cays UE1 of the vector mesons, without any re- 
adjustments of the model parameters. Hadronic 
processes involving four (or more) external par- 
ticles can also be described in this framework, 
though one has to go beyond impulse approxi- 
mation. Our results for ir-ir scattering agree 
with predictions from chiral symmetry, and this 
method has been applied to the anomalous "firirir 
coupling •">() as well. This approach unambigu- 
ously incorporates bound state effects in pro- 
cesses that receive contributions from off-shell in- 
termediate mesons, such as a and p mesons in 
case of 7T-7T scattering 49 , and VMD effects in 
the case of electromagnetic interactions. 

Within the same framework, one can consider 
baryons as bound states of a quark and a diquark, 
the latter in a color anti-triplet configuration. It 
has been shown |SJ that one can obtain a rea- 
sonable description of the ground state octet and 
decouplet baryons using both scalar and axial- 
vector diquarks. There has also been significant 
progress in understanding the nucleoli form fac- 
tors using this approach 52. 
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